SINP-APC- 11/01 



The Real Gauge Singlet Scalar Extension of Standard 
Model: A Possible Candidate of Cold Dark Matter 



Anirban Biswas §, Debasish Majumdar § 

Astroparticle Physics and Cosmology Division, 
Saha Institute of Nuclear Physics, 
1/AF Bidhannagar, Kolkata 700064, India 



ABSTRACT 

We consider a simplest extension of Standard Model in which a real SM gauge singlet scalar 
with an additional discrete symmetry Z2 is introduced to SM. This additional scalar can be a 
viable candidate of cold dark matter since the stability of S is achieved by the application of 
Z2 symmetry on S. Considering S* as a possible candidate of cold dark matter we have solved 
Boltzmann's equation to find the freeze out temperature and relic density of S for Higgs mass 120 
GeV in the scalar mass range 5 GeV to 1 TeV. As HHSS coupling 62 appearing in Lagrangian 
depends upon the value of scalar mass ms and Higgs mass nih, we have constrained the nis — 62 
parameter space by using the WMAP limit on the relic density of dark matter in the universe and 
the results of recent ongoing dark matter direct search experiments namely CDMS-II, CoGeNT, 
DAMA, EDELWEISS-II, XENON-10, XENON-100. From such analysis we find two distinct 
mass regions (a lower and higher mass domain) for such a dark matter candidate that satisfy 
both the WMAP limit and the experimental results considered here. We have estimated the 
possible differential direct detection rates and annual variation of total detection rates for this 
scalar dark matter candidate 5* for two detector materials namely Ge, Xe. Finally we have 
calculated the 7— ray flux from the galactic centre due to annihilation of two 130 GeV scalar 
dark matter into two monoenergetic 7— rays. 
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1 Introduction 



In recent years, one of the most important areas of modern cosmology is to investigate the 
existence and nature of dark matter in the universe. The observations by Wilkinson Microwave 
Anisotropy Probe (WMAP) |JJ for studying the fluctuations in cosmic microwave background 
radiation reveal that the universe consists of 27% matter and the rest 73% is an unknown energy 
known as Dark Energy. Out of this 27%, only about 4% accounts for the ordinary matter like 
leptons and baryons, gas, stars and galaxies etc. The rest about 23% matter is completely 
unknown. Moreover there are several cosmological observations like rotation curves of spiral 
galaxies, the gravitational micro-lensing, observations on Virgo [2] and Coma clusters [3j, bullet 
clusters [1], etc. which provide indications of the existence of huge amount of non-luminous 
matter or dark matter (DM) in the universe. 

Nature and identity of the constituents of dark matter are mostly unknown. However, evi- 
dences suggest that the dark matter candidates are mostly stable, non-baryonic, massive, non- 
relativistic particles having negligible or very weak interactions with other particles. These types 
of dark matter are often termed as cold dark matter (CDM) or weakly interacting massive parti- 
cles (WIMP). In the early universe, these particles would have been present in large numbers in 
thermal equilibrium. As the universe expands and cools down their density decreases resulting in 
decrease in their interaction/annihilation rates. When the expansion rate of the universe becomes 
larger than the annihilation rate of the WIMPs, they get decoupled from the universe. Thus 
they "freeze out" from the other contents of the universe and remain as relics. The temperature 
at which this phenomenon occurs is known as "freeze out" temperature and its density is called 
"relic density". After freeze out, the relic density of WIMP is only affected by the expansion of 
the universe. Since Standard Model (SM) of particle physics cannot provide any viable candidate 
for cold dark matter, one has to consider theories beyond SM in order to explain the dark matter 
candidates (namely WIMP). 

In this paper we have considered the simplest possible renormalisable extension of SM by 
adding a real gauge singlet scalar S. We impose a discrete symmetry Z2 on S and due to 
this symmetry the additional scalar S is stable and can be a viable candidate for cold dark 
matter. This model was first proposed by V. Silveira and A. Zee [5]. Thereafter a number of 
authors have explored its phenomenology [6]. The relevance of the scalar singlet as a plausible 
candidate for dark matter is very elaborately described in Ref. [7] (and references therein). 
Investigating the relic density of a scalar dark matter by constraining the unknown parameters 
from direct detection experiments are addressed by previous authors. In Ref. [S] , the relic density 
is investigated for scalar singlet by constraining dark matter mass and direct detection rates from 
DAMA [9] results. Similar analysis including the CoGeNT ^TO] results and CDMS II ^11] results 
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are also addressed in Ref. [12]. The analysis of scalar singlet dark matter scenario for XENON 
100 [13] direct detection experiment results are also given in this reference. The scalar singlet 
dark matter with CoGeNT results are also discussed by Fitzpatrick et al [H]. The interpretation 
of Fermi-Lat results with scalar singlet dark matter is discussed in Ref. [TB] . 

In the present work we estimate the freeze out temperature and relic density of the dark 
matter candidate S by solving Boltzmann's equation. Then we constrain the parameter 62 § 
by using WMAP limit on relic density of dark matter and the results of recent dark matter 
direct detection experiments like CDMS-II [IB HI], XENON-10 [IB], XENON-100 [IB], CoGeNT 
[ini US], EDELWEISS-II [2U3 and DAMA [S]. In CDMS and CoGeNT experiments the target 
material is Ge and in XENON experiments the target materials is Xe. 

The constrained parameters thus obtained are then used to calculate the differential direct 
detection rates and the annual variation of total detection rates of the scalar dark matter can- 
didate S for two detector materials namely Ge, Xe. Therefore we have calculated the 7— ray 
flux due to 130 GeV scalar dark matter for the annihilation channel SS — )• 77 from the galactic 
centre. 

The paper is organised as follows. In Section 2 we give a brief description of the scalar 
singlet model. Section 3 describes the formalism for computing relic abundance of a particular 
particle candidate. The results of the relic density calculations are given in Section 4. The model 
parameter 62 is constrained using the WMAP relic density data and the results obtained from 
various dark matter direct detection experiments. This is described in Section 5. The formalism 
for the calculation of direct detection rates and the annual variations of these rates is described 
in Section 6. With the constrained model parameter, 62 as obtained in Section 5, the direct 
detection rates and their annual variations of total detection rates are calculated for this scalar 
dark matter candidate for some reference detector materials namely Ge, Xe. The calculational 
procedure and the results are described in Section 7. In section 8 we have calculated the 7— ray 
flux from galactic centre due to annihilation of dark matter present in the galactic halo. Finally 
in Section 9, we give a summary and conclusion. 



2 The Model 

In the present work we consider a simplest extension of Standard Model where a real singlet 
scalar is added to the scalar sector of SM and explore the possibility that it can be a candidate 
for cold dark matter. The most general form of the potential appearing in the Lagrangian density 

^only parameter in this model which appears in both the expressions of scattering and annihilation cross 
section of S and which depends on the masses of scalar S and Higgs h 
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for such a scalar fields is 



V{H, S) = ^H^H + \{H^Hf + |ift^5 + ^Ah^hS' + j^^^ S +^-^S' + + ^-±S' 



and the Lagrangian of this model is given by 



(2) 



Where £sm is the Standard Model (SM) Lagrangian, H is the SM Higgs doublet and S is the 
real gauge (SU(2)lxU(1)y) singlet scalar. The stability of S will be achieved by imposing a Z2 
symmetry {S —S, C ^ C) over S. Therefore, under this symmetry the coefficients of odd 
powers of S are zero (/cs and 5i in Eq. ([2])). After spontaneous symmetry breaking masses of the 
Scalar field S and physical Higgs h are 

ml = fc2 + (52VV2, (3) 



ml = -m^ = AVV2 , (4) 

V is the VEV of Higgs (V = 246 GeV). In the present work we have taken the mass of the 
scalar particle S in the range 5 GeV - 1 TeV. Depending on its mass the dark matter candidate 
S annihilates into fermion pairs, gauge boson pairs and Higgs pairs. 



3 Formalism for Calculation of Relic Abundance 

In order to calculate the relic abundance of the scalar particle S we have solved numerically the 
Boltzmann's equation which is given by 

^ + 3Hn = -M(n2-ny , (5) 

where n is the number density of the scalar particle S and Ueq is the value of n when S was in 
equilibrium (when temperature T > Tf,Tf being the freeze out temperature of S*), H denotes 
the Hubble parameter, (av) is the thermal average of the product of annihilation cross section 
and the relative velocity of the two annihilating particles (in this case the scalar singlet S). It is 
useful to define two dimensionless quantities, Y = n/s [21] and x = m/T. Where s is the total 
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entropy density of the universe, T being the photon temperature. From the standard Friedmann- 



Robertson- Walker cosmology, the Hubble parameter H = y ^TiGp and G is the gravitational 
constant. The total energy density (p) and the total entropy density (s) of the universe are given 
by [21] 

P = 9eff{T)^T' (6) 

and s = Kff{T)—T' . (7) 

In Eqs.([n]) and ([7j) Qeff, ^e// are the effective degrees of freedom for the energy and entropy 
densities. Substituting Eqs. ([6]), ([7]) and the expression of H into Eq. ([5]), we arrive at the 
equation for the evolution of Y as 

where gl^"^ is defined as [21] 

Yeq is the value of Y when n = rzgq. The expression for Yeg is given by [21j 

Air'^Kjfim/x) ^ ' 

where g is the number of internal degrees of freedom of the species under consideration (here g 
= 1), m is the mass and Kn{x) is the modified bessel function of order n. From Eqs. ([8]) and 
ffTU]) . we obtain 

/45 \"'/'45^ K2{x) , ,,,, Ki{x) ld\nh,{T) 

In the above equation, hc{T) is the contribution to heff{T) from all species which are coupled to 
the universe at temperature T. Eq. lfTT]) above is solved numerically in a self consistent manner 
in order to obtain the value of Xf (and hence the freeze out temperature Tf (= m/xf)). In the 
present case we have taken the value of 6 to be 1.5 [21J. Integrating Eq.(|H]) from x = Xq = m/To 
to X = Xf = m/Tf, where Tq is the present photon temperature which is of the order of 10"^"^ 
GeV (~ 0), we obtain Yq (value of F at T = Tq). Knowing Yq we can compute the relic density 
of the dark matter candidate (here S) from the relation [21], 

771 

nh^ = 2.755 X 10^- — Yq . (12) 

GeV 
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In the above VL = pjpc (p being the dark matter density and pc is the critical density of the 
universe) and h = ^oo Km stf-iMpc-^ • Feynman diagrams in Fig. [1] represent the possible annihi- 
lation channels of S. The expressions for annihilation cross section {av) for the processes like 
SS — 7- //, W^W~ , ZZ, hh are given in refs.[221 123] In this work we consider Higgs mass value 




r \ ✓ V 



Figure 1: Lowest order Feynman diagrams of two S annihilate into a pair of fermion and anti- 
fermion , W^W~ , ZZ and Higgs. 

ruh = 120 GeV. The variations of annihilation cross sections with scalar mass ms are shown in 
Fig. [2] (Left Panel) for different values of 62- 



4 Calculational Procedure and Results 

The relic density for scalar dark matter is obtained after an elaborate computation. We first 
calculate the freeze out temperature Tf for scalar dark matter with different values of coupling 
constant 62 and mass ms- For this purpose we have solved Eq. f llip numerically. The values of 
the quantities gl^"^, glj^j- and hgff for different T required for solving Eq.f jTTj) . are obtained from 
the figures (for the QCD phase transition temperature of 150 MeV) given in Refs. [211 [2l]. Iii the 
Fig |5] (Right Panel), representative plots are the variations of Tf in the scalar dark matter mass 
range 5 GeV to 1 TeV for different values of 62 (^2 = 0.05, 0.1, 0.7) and nih = 120 GeV with the 
topmost plot is for the lowest value of 62 considered and the plots below are for other considered 
values of 62 in the increasing order. In general, the freeze out temperature Tf is approximately 
given by Tf ~ ms/ 20. The plots for Tf = 7715/20 are also shown in Fig. [2] (Right Panel) (black 
dashed lines) for reference. The sudden dip in the values for Tf in Fig. [2] (Right Panel) around 
ms = 60 GeV can be understood from the expression of {crv) fj (given in refs. [221 E3]). At 
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Mass of Scalar Dark Matter (Gev) Mass of Scalar Dark Matter (GeV) 

Figure 2: Left Panel : Variation of product of annihilation cross section and relative velocity 
{av) with the mass of scalar dark matter S for 82 = 0.7, 0.3, 0.1, 0.05, Right Panel : Variation 
of freeze out temperature Tf with the mass for different values of 62 = 0.7, 0.1, 0.05 



this point there is a sudden rise in (av) (Fig. [2] (Left Panel)) and it is due to Higgs propagator 
appearing in the annihilation process {SS — t- //). Using the values of freeze out temperatures 
(calculated from Eq. flTT]) ) in Eq. f[T2]) the relic densities of the scalar dark matter for different 
scalar dark matter masses and different values of 62 are computed. The results are shown in 
Fig. [3l In Fig. [3] the two parallel lines denote the WMAP limits on relic density of dark matter 
(0.099 < Qh"^ < 0.123). The different plots in Fig. [3] correspond to different values of 62 namely 
62 = 0.05, 0.1, 0.3, 0.7 respectively with the topmost one is for smallest value of 62 considered 
and the successive lower plots are for the other considered values of 62 in increasing order. We 
have seen from Fig. [2] (Left Panel) that initially the annihilation cross section of S increases 
with rris; then at m5 ~ iTih/^, (av) rises rapidly and after which it decreases with the increase of 
rris- Again for 7715 ~ 81 GeV, (av) suddenly increases upto nearly 2 orders of magnitudes from 
its value at 777,5 ~ 80 GeV and it is due to the fact that for ms > 80.4 GeV the annihilation 
channel SS — )■ W^W~ becomes kinematically possible. Thereafter (erf) starts decreasing with 
the increase of ms- Since relic density is inversely proportional to {av) 0, the variation of relic 
density of dark matter particle S with is just opposite to the variation (av) with ms- This 
feature is reflected Fig. [31 Also since {av) is directly proportional to 5| and its higher powers, 
higher the value of 82 lower is the value of relic density (Fig. [3]). 

Physically we can say that {av) is directly proportional to probability of that process. So for higher {av) 
the probability of pair annihilation of S is high and hence density is low 
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S2 = 0.7 
WMAP Lower Limit 
WIVIAP Upper Limit 



10^ 10^ 10^ 

Mass of Scalar Darl< Matter (GeV) 

Figure 3: Variation of relic density Qh'^ with the mass of scalar dark matter for Higgs mass 120 
GeV 

5 Constraining the model parameter 62 

The model parameter 62 is a very important parameter for this present model because it appears 
in both the expressions of annihilation and scattering cross section of scalar dark matter S. 
The spin independent scattering cross section for scalar dark matter S is given later in Eq. 
(fT3|) . In this section we have constrained the parameter space {ms — 62) by using WMAP limit 
on relic density of dark matter and the results of recent experiments like CoGeNT, DAMA, 
CDMS-II, XENON-10, XENON-100 and EDELWEISS-II. Similar to the previous discussions, 
here also we perform the calculations for Higgs mass rrih = 120 GeV with ms in the range 5 
GeV < nis < 1 TeV. The results obtained are shown in Fig. HI We first use the WMAP limit 
(0.099 < Qh"^ < 0.123) on relic density of dark matter and using that limit we get the allowed 
values of for each value of 62 (from Fig. [3]). These results are shown in Fig. H] using turquoise 
coloured contour. Thereafter we estimate the allowed values of 62 and ms using the mass - cross 
section hmits given by the experiments like CDMS-II, DAMA, CoGeNT, XENON-10, XENON- 
100, EDELWEISS-II and Eq. 0131) . In Fig. |3]the magenta coloured contour represent the allowed 
regions of scalar dark matter S obtained from CoGeNT data (2010). The overlap regions between 
these contours (magenta and Turquoise) are therefore satisfied by both WMAP and CoGeNT 
(2010) results. From the overlap region (Fig. H]) the range of ms (in GeV) is found to be 
7.7 < ms < 11.15 and the corresponding range of coupling 62 is obtained as 0.7 < 62 < 0.95. 
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Mass of Scalar Dark Matter (Gev) 

Figure 4: Constraining the parameter space {ms — 62) of scalar dark matter for Higgs mass 
120 GeV (upper panel) by using WMAP limit and recent experimental results of CDMS-II, 
XENON-10, XENON-100, CoGeNT, DAMA, EDELWEISS-II. 

These ranges therefore satisfy both WMAP and CoGeNT (2010) limits. The dark blue coloured 
contours indicate (Fig. H]) new bounds from CoGeNT data (2011) flQlj. The common region 
between CoGeNT (2011) and WMAP lies in the range 8.0 < mg < 8.53 GeV, 0.8 < ^2 < 0.9. 
This common intersection region is also well supported by CoGeNT (2010) and CDMS-II (2010) 
data [25] (purple dashed line). 

Similar ms — 62 contours obtained from the DAMA experiment results (with channeling) 
are shown as maroon contours in Fig. |H One sees that the small overlap regions between 
the two contours (turquoise and maroon) are restricted by the scalar mass (in GeV) range 
14.8 < ms < 15.9. The corresponding values of 62 are found around 0.6. We remark in the 
passing that we have checked for other allowed regions in "dark matter mass — o"nucieon" plane 
given by DAMA experiment but we have not obtained any overlap region such as described above. 
One of such regions (maroon coloured contour) is also shown in Fig. |H The olive dashed line in 
Fig. m represent upper bounds that we have obtained from low energy analysis of the CDMS-II 
Germanium data But it has no intersection with WMAP satisfied region (turquoise coloured 
contour) . 

But unlike CoGeNT and DAMA, other experiments like XENON-10, XENON-100, CDMS-II, 
EDELWEISS-II do not provide a bounded allowed region in ms - ctn plane (un is the scattering 
cross section of dark matter and nucleon). Instead they provide upper bounds of scattering cross 
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section for a particular mass of dark matter. Consequently we also obtain upper bounds of 62 
for a specific mass of S for those experiments. These results are also shown in Fig. |H In Fig. H] 
green dashed line represents the upper bound of 82 for XENON-100. CDMS-II results (CDMS 
2009 Ge and CDMS Soudan (All) @) are shown by black and blue dashed lines. XENON-10 
and EDELWEISS-II results are represented by red dashed, orange solid line respectively. The 
WMAP results (turquoise plot in Fig. H]) intersect with the upper bounds obtained from XENON- 
10 results in ms — 82 plane are found to be at the values of ms = 53.5, 62 = 0.12, ms = 67.4, 
62 = 0.16 and ms = 80.2, 62 = 0.21. Therefore the turquoise region below {ms = 53.5, 62 = 0.12), 
(ms = 67 A, 62 = 0.16) and {ms = 80.2, 62 = 0.21) is obeyed by both WMAP and XENON-10. 
Similarly the regions satisfied by both CDMS-II results (CDMS 2009 Ge, CDMS Soudan (All)) 
and WMAP are represented by the turquoise colour below intersection points {ms = 52.0, 62= 
0.15), {ms = 67.9, 82 = 0.17), {ms = 80.3, 82 = 0.20) and {ms = 53.9, 82= 0.12), {ms = 66.5, 62= 
0.14), {ms = 80.4,^2= 0.16) respectively as shown in the Fig. HI Also the overlap regions of 
WMAP, XENON-100 and WMAP, EDELWEISS-II are below the following intersection points, 
can be read out from the Fig. |4]as {ms = 57.6, 62 = 0.05), {ms = 62.5, 62 = 0.05), {ms = 80.5, 
62 = 0.07) and {ms = 54.1, 62 = 0.11), {ms = 66.7, 82 = 0.14), {ms = 80.4, 62 = 0.17). We have 
also found that for XENON-100 there is another intersection point with WMAP in the lower 
mass region around ms ~ 6.0 GeV, 62 ~ 1.25. 

From the above analyses it appears that there are two distinct regions in the ms — 62 plane 
for scalar dark matter 5* which are allowed regions for both WMAP and recent experiments. The 
regions can be classified as follows. 

• A lower mass region where we have found 3 mass ranges for scalar dark matter S. These 
ranges are given by ~ 6 GeV (^2 ~ 1-25), 7.7 GeV < ms < 11.15 GeV (0.7 < ^2 < 0.95) 
and 14.8 GeV < ms < 15.9 GeV (^2 ~ 0.6). The corresponding ranges for coupling 62 
which we have found are given within brackets. This lower mass domain is supported by 
WMAP and various ongoing dark matter direct detection experiments. In this case ms ~ 6 
GeV is supported by WMAP and XENON-100. Second and third mass ranges are obeyed 
by WMAP, CoGeNT (2010) data and WMAP, DAMA (with channehng) data respectively 
(Fig. SD. 

But if we use more recent data of CoGeNT (CoGeNT (2011) data) then the second mass 
range of scalar dark matter S gets reduced to 8.0 GeV < ms < 8.53 GeV. The ranges for 
coupling 62 also reduced to 0.8 < ^2 < 0.9. It is also seen from Fig. H] that is region is 
supported by CDMS-II (2010) bounds. Other mass ranges are remain unchanged. 

^Which are the results obtained by the CDMS-II collaboration from the combined analysis of full data set of 
Soudan. 
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• A higher mass region with the scalar dark matter mass range ~ 52.5GeV < ms < ~ 1000 
GeV, with the range of 62 found as 0.02 < ^2 < 0.4 (Fig. H]). This mass region is satisfied by 
the allowed domains of WMAP, CDMS-II, EDELWEISS-II, XENON-10 and XENON-100. 

Here we make some comments about the region of the parameter space {ms vs ^2) that is 
not satisfied by the results of direct detection experiments we have considered. In this region (15 
GeV< ms < 52.5 GeV) the values of 62 required to obtain current relic density (within WMAP 
limit) are such that the scattering cross sections obtained using these values for different dark 
matter masses do not satisfy the experimental limits given by different exclusion plots. 



6 Formalism for the Calculations of Direct Detection Rates 

In this section we estimate the differential direct detection rates and their annual variations for 
scalar dark matter S. For this purpose we have chosen '^^Ge and ^^^Xe as a detector materials. 
The direct detection of dark matter by a terrestrial detector uses the principle of elastic scattering 
of dark matter particles off the detector nuclei and the energy of the recoil nucleus is measured. 
It is very difficult to measure the low recoil energy of nuclei accurately and hence a very low 
threshold and low background detector is required. In Fig. [5] we show the Feynman diagram for 
such elastic scattering process of scalar singlet 5* through Higgs exchange. The scalar singlet S 
- nucleon elastic scattering {SN — > SN) cross section is given by 

scalar SyjA^ f m^\ 



Air \Ms^Mh' 

where, mr{N, S) = M^Ms / {Mn + Ms) is the reduced mass, A is coupling between Higgs and 
nucleon A^ and its value is ~ 340 MeV/V^ [26], with V being the VEV of Higgs boson. The scalar 
singlet - nucleus elastic scattering cross section is given by 



scalar ^ A^m^ (uUclcUS, S) scalar 
C^nucleus m2 (nUclcOU, ^) ' ^ ' 

In the above A is the mass number of the nucleus. The differential detection rate of dark matter 
per unit detector mass is given by 



J p ^scalar „ 

an _ Onucleusr5'^2^^^j 



(IEr Avem^sra. 

X 

TiT2T^ 



erf — erf 



(15) 
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Figure 5: Feynman diagram for the elastic scattering between 5* and nucleon N via Higgs ex- 
change. 



where 



T nucleusr t-i t-i2 / 7-1 

1 = Z ^) ±2 = l-C/, 



T- 
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eri — eri 



(16) 



Where Ve is the velocity of earth with respect to galactic frame of reference. And its expression 
is given by [27] , 

Ve = t;© + forbCOS7COS ( I . (17) 

In the above expression t denotes any time of the year, T = 1 year is the time period of earth's 
motion around the sun, Vqa = 30 Km/sec is earth's orbital speed and 7 ^ 60° is the angle 
subtended by the ecliptic at the galactic plane. The solar velocity Vq is given by 

Vq = Vq + Vpec , (18) 

where vq is the circular speed of sun around the galactic centre taken to be 220 Km/sec and Vpec 
is the peculiar velocity with Vp^c = 12 Km/sec. The periodicity in Eq. (1171) causes an annual 
modulation of the event rates of dark matter in a terrestrial detector which serve as a definite 
signal of dark matter detection. In the Eq. US] F (£'/?) is the nuclear form factor given by [28], 
ps is the dark matter density in the solar neighbourhood, equal to 0.3 GeV/cm^ for the rest of 
our calculations in this section . Vmin denotes the minimum velocity of dark matter required to 
produce a recoil energy E^. The expression of Vmin is given by. 



(19) 
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The measured response of the detector by the scattering of dark matter off detector nucleus 
is in fact a fraction of the actual recoil energy. Thus, the actual recoil energy Eji is quenched by 
a factor qx (different for different nucleus X) and we should express differential rate in Eq. f|T5|) 
in terms of E = QxEr. Thus the differential detection rate (events/Kg/Day/keV) in terms of 
the observed recoil energy E for a monoatomic detector like Xe can be expressed as 

(^) = / -77?(^^)^ • (20) 



The total detection rate of dark matter is obtained by integrating Eq.( lT^ as 

^ = / ^'^^^ ' ^21) 
JEt "-Er 

where Et is the threshold energy for a given detector material. 



7 Direct Detection Rates for Scalar Dark Matter 

In the present work, computations of direct detection rates are performed with = 120 GeV, 
AE = 0.5 keV and at a time t = to. As discussed earlier we have computed the direct detection 
rates and their annual variations for each of the detector materials namely ™Ge, ^^^Xe. The 
quenching factors for ^^Ge = 0.25 ^^^Xe = 0.8 [29], The differential detection rates and their 
annual variations can now be computed using Eqs. ffT^ - ([21]) • 




Figure 6: Variation of differential detection rates AR/AE of scalar dark matter 5* with observed 
recoil energy E for monoatomic detectors Xe (left panel), Ge (right panel) 

The variation of differential detection rates of scalar dark matter S with observed recoil 
energy E for mono atomic targets like Xe, Ge are shown in Fig. [61 In the left panel of Fig. [6l 
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we show the estimates of differential detection rates for different values of observed recoil energy 
E, with Xe as target material for ms = 55 GeV (red solid line) and ms = 65 GeV (green solid 
line). For both the cases the value of the coupling constant 62 is taken to be 0.1 (in agreement 
with the higher mass region described in section 5). Left panel of Fig. IH] shows that although the 
two plots corresponding to two scalar masses are distinguishable at lower recoil energies (< 11 
GeV), at higher recoil energies they tend to coincide. In the right panel of Fig. [6] we show the 
direct detection rates results for the case of Ge. In this case, calculations are performed for two 
sets of ms — 62 values namely {ms = 10 GeV, 62 = 0.8), (m^ = 8 GeV, 62 = 0.9). These values 
are chosen from the allowed lower mass domain discussed in section 5 (Fig. H]). It is seen from 
the right panel of Fig. [6] that the rates for the set {ms = 8 GeV, 62 = 0.9, represented by green 
solid line), falls off faster than those for the set {ms = 10 GeV, 62 = 0.8, represented by red 
solid line). The nature of AR/AE for the case of Xe (left panel of Fig. |6]) can be explained 




Figure 7: Left panel - Variation of T3 with observed recoil energy E for Xe with ms = 55 GeV 
(green solid line), 65 GeV (red solid line). Right panel - Same for Ge with ms = 10 GeV (red 
solid line), 8 GeV (green solid line) 

by examining the variation of T3 (Eq. [TB]) with E. This is shown in left panel of Fig. [71 For 
low values of E and for high mass range of S {ms > 55 GeV), Vmin (Eq. [I9])<^ Ve and hence 
T3 is effectively independent of E. Therefore as E increases, the values of T3 for ms = 65 GeV 
becomes larger than those for ms = 55 GeV. Also cr^^'^^gyg is inversely proportional to ms and 
Ti is directly proportional to cr^'^'^^^^g. Consequently Ti is inversely proportional to ms- Now the 
variation of AR/AE with E is due to the combined effects of both Ti and T3 (Eq. f ll5p ). This 
explains the nature of the plots for Xe in Fig. |6l In the case of Ge however, Ti is nearly the 
same for both the masses considered. Consequently the variations of AR/AE with E (Fig. |6] 
(right panel)) is dominated only by the nature of variations of T3 with E. This variations are 
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shown in right panel of Fig. [7] which explains the nature of variations of AR/AE with E for Ge. 

The annual variations of total detection rates of WIMP is a crucial evidence for dark matter. 
This variation is caused by the periodic motion of earth around the sun in which the directionality 
of earth's motion changes over the year. Since the solar system moves towards the direction of 
Cygnus constellation, earth experiences a WIMP wind apparently coming from the direction 
of Cygnus. Due to the periodic motion of earth, the relative speed between earth and WIMP 
changes over the year. It becomes maximum when both the velocities of solar system and earth 
are in the same direction (on 2^'^ June) in which case the earth encounters maximum WIMP 
flux. The WIMP flux encountered by the earth is minimum when velocities of earth and sun are 
in opposite direction. Consequently, maximum events are expected on 2"'^ June of every year. 




50 100 150 200 250 300 350 400 50 100 150 200 250 300 350 400 



Figure 8: Left panel - Annual variation of total detection rates of scalar dark matter for Xe 
(mono atomic target) with ms = 55 GeV, 62 = 0.1. Right panel - Same for Ge with ms = 10 
GeV, 62 = 0.8 

In this work we compute the total detection rates at each day of a year, for the same detector 
materials namely Xe, Ge considering the scalar singlet as dark matter candidate. The results are 
then plotted with the days of year which show the annual variation of total detection rates. The 
calculations for Xe are performed for the set {ms = 55 GeV, 62 = 0.1) whereas for Ge, the set 
{ms = 10 GeV, 62 = 0.8). The results for Xe and Ge are shown in left and right panels of Fig. 
E] respectively. All the plots in Figs. Elshow that the maximum expected events are at t = 153 
(day) (corresponds on 2"^^ June). 
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8 Indirect Detection of Scalar Dark Matter 



Another promising method for the detection of dark matter (WIMPs) is the observation of 
annihilation products of dark matter present in the galactic halo. In this section we will consider 
7-rays coming from the dark matter annihilation in the galactic centre (GC). 

Recently it has been reported that there is a 4. 60" (3.3cr) [30], |3T] local (global) evidence of a 
monochromatic gamma-ray line with an energy ^ 130 GeV by the publicly available data [32] 
of Fermi Large Area Telescope (Fermi-LAT). This signal comes from two regions one of which 
is nearly at the centre of our galaxy (—1", —0.7°), hereinafter referred to as the "central region" 
and another is located at (—10°, 0°), called the "west region". Both regions are extended with in 
a circle of radius of 3°. It is suggested that this excess of gamma ray signal from galactic centre 
(GC) is not associated with the Fermi bubbles [21] and may result from dark matter annihilation 
into two monochromatic gamma-rays. 

We have calculated the 7-rays flux due to 130 GeV scalar dark matter annihilation in the 
"central region" of our Milky way galaxy. The Feynman diagram for the process SS — > 77 is 
shown in Fig. [9l The expression of 7— ray flux due to dark matter annihilation in galactic halo 




>---*- — 

✓ 

/s 



Figure 9: Feynman diagram for the process SS — 77 

is given by [33j, 

-rePiJ , (22) 



dE^ 8tt rrii dE^ 



where 



ds .fp{ry^ 



and 



J = I db I dl — cos6 ^ ) (23) 
'i.o.s \ Pe 



^ = 26{E - E,) . (24) 
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Mass of Scalar 
Dark Matter 

GeV 


Coupling 82 


Flux using 
Einasto Profile 

GeV cm~^ s^-*^ sr^-*^ 


Flux using 
NFW Profile 

GeV cm^^ s^^ sr^^ 


Flux using 
Isothermal Profile 

GeV cm~^ sr~^ 


130 


0.06 


1.971 X 10-^ 


9.801x10-8 


4.048 X 10-9 



Table 1: 7— ray flux obtained from the annihilation channel SS — )■ 77, for three different dark 
matter halo profiles. 

In the above, / and h are the galactic longitude and latitude respectively. We have performed /, 
b integration (in Eq. [23|l over the "central region" of the our galaxy and the s integration (in 
Eq. [2^ along the line of sight (l.o.s). Relation between r and s is given by 

r = (s^ + Tq — 2srQ cos/ cos 6) 2 , (25) 

where = 8.5 Kpc, the distance of the sun from the galactic centre and = 0.4 GeV/cm'^ is the 
dark matter halo density at the position of the solar system. The expression of the annihilation 
cross section {(tv) ^^^^^ (in Eq. 22) for the process shown in Fig. [Hlis given in Ref. [31]- In this 
calculation we have taken three different dark matter halo profiles (available in literature) namely 
the Einasto profile [35], the NFW profile [36] and the Isothermal profile [37]. These halo profiles 
give the functional dependence of p(r) with r. In the present calculation we have considered the 
value of the coupling 82 = 0.06, which is allowed by WMAP and all recent ongoing dark matter 
direct detection experiments that have been considered in this work (Fig. IU section 5). We have 
calculated the 7— ray fluxes for all the halo three profiles considered above and the results are 
shown in Table 1. The annihilation cross section (av) _^ is calculated to be 7.13 x 10"^^ cm^/s 

\ / SS — / 

for 62 = 0.06. From Ref. [31], one sees that the 7-ray flux obtained from "central region" of our 
galaxy is in the range 4.0 x 10"^ to 7.5 x 10"^ (in GeV cm-^s'^sr"^) (95% CL) with best fit value, 
5.6 X 10"^ GeV cm-^s-^sr"^. From Table 1 we see that in order to compare our results to those 
in Ref. [31] the annihilation cross section for the channel 5*5* — 77 in the present calculation 
must be enhanced by a factor of ~ 3.0 x 10^ (for the Einasto profile) and ~ 5.7 x 10^ (for the 
NFW profile) with respect to best fit value. As a result we have to increase coupling 62 from 
0.06 to 1.03 (for the Einasto profile) and 1.43(for the NFW profile) respectively. From Fig. Hjit 
is seen that such a high value of 62 is not satisfied by either any direct detection experiments we 
have considered in this work or by the WMAP limits. Therefore we conclude that a 130 GeV 
dark matter in the present model can not explain the Fermi-LAT observed 130 GeV 7— ray line 
in the direction of the galactic centre unless the process is boosted (by introducing a boost factor 
[38]) either by astrophysical justifications and/or by other particle physics methods. 
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9 Summary and Conclusion 



In the present work we consider a simplest extension of SM by introducing a real gauge singlet 
(singlet under SU(2)lxU(1)y) scalar S to SM which can only interact with SM particles via 
Higgs. For the stability of S, Z2 symmetry is imposed in the theory. Thus S can be a viable 
candidate for cold dark matter. The scalar mass ms and the coupling are the two parameters 
in the theory. We have calculated the freeze out temperature and relic density of this scalar 
dark matter candidate S by solving Boltzmann's equation and have constrained the ms — 62 
parameter space by using WMAP limit on relic density of dark matter in the universe and the 
results of recent ongoing dark matter direct search experiments like CDMS-II, DAMA, CoGeNT, 
XENON-10, XENON-100. We find that if S" is a dark matter candidate then its mass appears to 
be constrained within two regions. One is a lower mass region where mg can vary from 6 GeV to 
16 GeV with 62 lies in the limit 0.7 < ^2 < 1.25 for ruh = 120 GeV. This region is supported by 
WMAP, CoGeNT and DAMA data. The other region is higher mass region with the ranges for 
ms (in GeV) and 62 found out to be 52.5 ^ ms ^ 1000, 0.02 < 62 < 0.4 for the same Higgs mass. 
This region is also supported by the limits given by WMAP, CDMS-II, XENON-10, XENON-100, 
EDELWEISS-II experiments. We have calculated the possible differential direct detection rates 
and annual variations of total detection rates for scalar dark matter S in case of two detector 
materials namely Ge, Xe. For all these target materials we have found that differential detection 
rates decrease rapidly with the increase of observed recoil energy and they become vanishingly 
small for recoil energies beyond 10 keV for Ge with scalar mass = 10 GeV. Whereas for Xe, 
the rates become vanishingly small for recoil energies beyond 80 keV when = 55 GeV. We 
have also shown how the total rates vary over a year for these target materials. These annual 
variations of total detection rates, if found, will be one sure evidence for dark matter detection. 
Finally in the last section we have calculated the 7— ray flux for a 130 GeV scalar dark matter S 
and we have found that it is not possible to explain the Fermi-LAT observed excess 7— ray line 
by a 130 GeV scalar dark matter, unless a boost factor of order of 10^ is introduced with the 
annihilation cross section of SS — )■ 77 channel. 

Acknowledgments: A.B. thanks Debabrata Adak for some valuable discussions. 
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